Phonon-induced topological transitions and crossovers in Dirac materials 
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We show that electron-phonon interactions can alter the topological properties of Dirac insulators 
and semimetals, both at zero and nonzero temperature. Contrary to the common belief that in- 
creasing temperature always destabilizes topological phases, our results highlight instances in which 
phonons may lead to the appearance of topological surface states above a crossover temperature in 
a material that has a topologically trivial ground state. 
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Introduction. — Recent years have witnessed the un- 
expected discovery of novel topological phases in cer- 
tain electrically insulating [l[ and semimetallic [|| solids 
whose low-energy bulk excitations are Dirac fermions. 
Topological phases in these "Dirac materials" are char- 
acterized by nonzero integers (topological invariants), 
which manifest themselves through peculiar and robustly 
gaplcss states localized at the boundaries of the material. 
Partly enticed by a vision of transistors that would oper- 
ate by switching topological invariants on and off, there is 
keen interest in finding ways to induce topological phases 
in intrinsically non-topological materials. 

In ccntrosymmctric crystals, quantum phase transi- 
tions between topological and non-topological insulators 
occur when the sign of the Dirac fermions' mass gets in- 
verted at an odd number of time-reversal-invariant mo- 
menta (TRIM) in the bulk Brillouin zone. In equilib- 
rium, such "band inversions" can result from applying 
pressure 0| or dc electric fields from changing the 
compound stoichiometry or from strong alloying Q. 
Away from equilibrium, ac electromagnetic fields may 
effectively produce a band inversion [7(, and engineered 
coupling to dissipative baths presumably induces topo- 
logical phases in driven cold atomic systems The 
latter two cases illustrate how the topological invariant 
for an electron system may be altered after coupling it to 
a non-electronic environment, and motivate the question 
of whether analogous environment-induced effects could 
occur in equilibrium. 

In crystals that are in thermodynamic equilibrium, lat- 
tice vibrations form the most ubiquitous bath that cou- 
ples to electrons. Although electron-phonon coupling can 
hardly be engineered in a given material, its influence on 
electronic properties is strongly temperature-dependent. 
Even as a recent formal study suggests that topological 
invariants are not robust under dephasing Q, a growing 
body of work (TMl treats phonons as mere actors that 
do not alter the topology of the underlying electronic 
stage. The objective of our paper is to demonstrate that 
(i) phonons can modify topological properties of an elec- 
tronic structure, and (ii) temperature can be used as the 
"agent" that drives those changes. 

Electron-phonon interactions are known to often de- 



crease the bandgap of ordinary semiconductors mono- 
tonically as a function of temperature [i~6j . However, it 
has not been recognized that, at least in Dirac materi- 
als, phonon-induced renormalization of the bandgap can 
culminate in a band inversion both at zero and nonzero 
temperature. In insulators with non-topological ground 
states, a band inversion arising at nonzero temperature 
implies the emergence of peculiar gapless surface states 
above that temperature. This may be regarded as a 
"topological thermal crossover" , which has not been pre- 
viously discussed. 

Electron-phonon self-energy.- Consider electrons in a 
perfect periodic lattice with crystal Hamiltonian ftA°)(k), 
whose energy eigenvalues Ek a and eigenstates l^ka) = 
exp(ik-r)|ka) characterize the electronic structure. Here 
r is the position vector, k is the crystal momentum, a = 
1,...,JV is the band index and |ka) is a A-component 
cigenspinor. The change in the electronic structure due 
to lattice vibrations is approximately captured by the 
electron-phonon self-energy [l?! [H[ 
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where w„ = (2n + 1)ttT (n € Z) is the Matsubara fre- 
quency, T is the temperature, £k Q = E^ a — cf, £f is 
the Fermi energy, w q is the phonon dispersion, n q = 
[exp(cj q /T) — 1] _1 is the phonon occupation number, 
jka = [exp(^ka/7 1 ) + l] _1 is the fermion occupation num- 
ber, and g q is the electron-phonon coupling. For con- 
creteness and simplicity we concentrate on the deforma- 
tion coupling to longitudinal acoustic phonons (l7j . so 
that (?q = hC 2 q/(2pa 3 c s ) and w q = c s q, where C is the 
deformation potential, a is the lattice constant, p is the 
atomic mass density and c s is the sound velocity. 

In presence of lattice vibrations, the electronic Green's 
function G obeys (iu; n , k) = (iu n - £ka) 6 aa > - 

T, aa r (iu> n , k), where 6 is Kronecker's delta. The off- 
diagonal matrix elements of the self-energy are gener- 
ally nonzero and lead to avoided band crossings ex- 
cept at certain (e.g. topologically protected) degener- 
acy points. Upon diagonalization, the Green's function 



2 



obeys G a ^(iuj ni k) 
the solutions of 



(iu r , 



e F 



Xa(«w„,k)) 6 aa >, and 



Et 



iTi, 



XaiE^a - £F + ir kQ ,k) 



(2) 



yield the energy dispersion E£ a and broadening |rk Q | of 
dressed quasiparticles at temperature T. 

Phonon-induced topological transitions in insulators. - 
A minimal lattice model that captures the generic low- 
energy physics of 3D Dirac insulators reads [19( 

ft<°>(k) = e (k)l 4 + d(k) • crr x + M(k)l 2 T z , (3) 

where 1„ is an n x n identity matrix, a 1 and t 1 are 
Pauli matrices in "spin" and "orbital" space (respec- 
tively), e = 27 (3 - J2i cos(fcia)), rf, = -2Asin(fcja), 
M(k) = m + 2t(3-X)iCos(fcio)), i £ (7, A,*) 

are band parameters and m is the Dirac mass. Anticipat- 
ing that \m\ < \t\ in typical Dirac materials, this model 
describes an insulator with an energy gap E g = 2m at 
k = 0: it is topological if mt < and non-topological 
if mi > 0. It contains two degenerate conduction 
bands (E^i = E^ 2 ) and two degenerate valence bands 
(£ k3 = £k4 = -Ski). Since |£ kce | > |m| for k ^ 
TRIM, phonon-induced changes in topological properties 
are realized via band inversions at k = 0. Thus we focus 
on Ti aa i (iu) n , 0), which can be recasted in matrix form as 



S(iw n , 0) = E («w„)l 4 + E z («w„)1 2 t z 



(4) 



with E (iu; n ) = (£n(m; n ,0) + S 33 (iw„, 0))/2 and 
E z (£w n ) = (Eii(iu> n ,0) - E 33 (ja; n ,0))/2. 

At T = 0, a topological quantum phase transition oc- 
curs when the renormalizcd Dirac mass, m* = m + £ z (0), 
crosses zero and changes sign • Due to the frequency- 
dependence of the self-energy, 2m* is in general not 
equal to the renormalized energy gap E* = E* — E*, 
where E* ~ m + S (£ ! * - e F ) + - e F ) and 

E* ~ -m + £ (£7* - £ f) - S Z (S* - £f)- However, in 
the toy model treated here, £ aQ ,/(w,k) ~ E aQ /(0,k) and 
thus £"* ~ 2m*. In other words, the sign change of m* 
occurs nearly simultaneously with the band inversion of 
the dressed bulk quasiparticle spectrum. 

One interesting aspect of phonons is that their influ- 
ence on the electronic structure is strongly temperature- 
dependent below the melting point of the solid, because 
the Debye temperature is of the order of room tempera- 
ture 2l|. Accordingly, m* is T-dependent and may in 
principle change sign as a function of temperature at 
some T c > 0. Provided that quasiparticle broadening 
is small compared to \m*(T)\, the topological proper- 
ties of G are identical to those of h^ ' with m = m*(T). 
Therefore, a sign change of m* (T) at T = T c implies the 
emergence (if m*(0)t > 0) or evanescence (if m*(0)t < 0) 
of topological surface states at T > T c . Angle- resolved 
photocmission may be the most direct probe to access 
this "thermal topological crossover" . 



Adopting the spherical approximation of Eq. (|3|) near 
k = and introducing a UV momentum cutoff k c = Tr/a, 
Eq. (j4j can be evaluated analytically in the limits T<w c 
and T>w c , where ui c = c s k c is the Debye frequency We 
obtain 



Eo, z =U 2 f , z (T), 



(5) 



where U = g cl /(aq) 1 ^ 2 is a constant in units of energy, 
and /o, z (T) are functions of electronic structure param- 
eters and of temperature, with 

f z (T « Wc ) ~ - y [1 - x 2 ln(l + l/x 2 )} /4 + 0(T 4 ) 
f z (T » u c ) ~ - (T/u B )y [1 - x tan" 1 (l/x)] , (6) 

x = {2/ir)\/{t 2 - 7 2 ) 1 / 2 and y = t/(t 2 - 1 2 ). Simi- 
larly, / (T) = -(7/*)/r(T). In the derivation of Eq. © 
we have neglected w q from the denominators of Eq. ([lj , 
which is a good approximation as confirmed numerically. 
This approximation leads to the cancellation of Fermi 
factors in Eq. ([I]) and thus highlights the fact that oj c 
governs the T-dependence of m* . 

From Eq. ([5]), a topological transition or crossover oc- 
curs at U — U C (T), where 



U C (T) = y/-m/f z (T). 



(7) 



Whether phonons favor topological or non-topological 
phases depends on the band parameters of the perfect 
crystal. If t 2 > 7 2 , sgn(/ z ) = — sgn(t) and phonons tend 
to drive the system into the topological insulator phase. 
If t 2 < 7 2 , phonons instead stabilize the non-topological 
insulating phase. Since |/ 2 (T)| is a monotonically de- 
creasing function of x, with f z — > as x — > 00, phonons 
are more likely to affect the band topology in materials 
with flatter bands (smaller A). 

Some reasonable parameter values are 
0.25 eV, 7 ~ 0.15 cV, A ~ 0.18eV, C ~ 



IK 



22J t ~ 
35 cV, p ~ 

7800 kg/m 01 , c s — 1.7km/s and a ~ lnm. Then, 
U ~ 70meV and J7 C (0) ~ 34i/rn]meVj. Accordingly, 
phonons may drive a topological quantum phase transi- 
tion in small-gap semiconductors (m < lOmeV). More- 
over, because U C (T 3> u> c ) <C U c (0), a material that is 
topologically trivial at T = may develop topological 
surface states at higher temperature. 

Aside from the Dirac mass, electron-phonon inter- 
actions renormalize the Fermi energy, e* F (T) ~ ep — 
£o(0) Nl- Though not explicitly, U c in Eq. Q depends 
on e* F implicitly because screening (ignored above) leads 
to # q 3 q /K(q), where «(q) ~ 1 + q\ F jq 2 and g TF 
is the Thomas-Fermi wavevector. Therefore, the obser- 
vation of phonon effects on topological properties may 
be contingent on having low bulk doping {qTF *C k c ). In 
view of this, electrically gated thin films are the preferred 
venue to tune the chemical potential inside the bulk gap 
and to observe transport and thermodynamic signatures 
of phonon-induced surface states. 
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FIG. 1: Dressed bulk quasiparticle dispersion along [001], 
obtained from Eqs. |T]), ^ and l[3)l. Parameter values are 
U ~ 57meV, 7 = 0.15 eV, A = 0.18 eV, t = 0.25eV, 
u) c — 8meV, ojq = w c /(1.24ir)o g and m = 20meV (where 
1.247r/a is the radius of the Debye sphere). Because electron- 
phonon interactions preserve lattice and time-reversal sym- 
metries, dressed bands are two-fold degenerate. Left: T — 
(topologically trivial insulator). Middle: T ~ 0.7 ui c (topo- 
logical crossover temperature). Right: T = 1.2 w c (inverted 
regime). The band broadening \Tit, a \ (not shown) grows with 
T, but remains < 1 meV throughout the three panels. 



The thermal topological crossover introduced above 
is meaningful only if the quasiparticle band broaden- 
ing in the vicinity of TRIM remains small compared to 
|m*|. The reduced phase space for scattering near k = 
helps keep the broadening small. For example, it fol- 
lows from inspection that there exists a finite window, 
1 S U 2 f z (T)/m < 2/(1 + 7/t), for which quasiparticle 
bands are inverted and yet sharp (Fo a = 0) at the same 
time. Here we have neglected the broadening due to e.g. 
impurities and bulk-surface scattering. The latter will be 
suppressed at k ~ insofar as \m*\ > (w c , T). 

Along similar lines, the topological surface states de- 
rived from G are hardly detectable unless \E*\ ~ 2\m*\ 3> 
T. Whether this condition is met or not depends on 
material parameters, because \m*\ oc T at high tem- 
perature (c.f. Eq. ©). In our toy model, using the 
parameter values enlisted above, \E*\ 3> T translates 
into uj c <C 25 meV, which can only be satisfied in soft 
Dirac materials with low Debye temperatures. This up- 
per bound on lo c becomes larger (and thus easier to fulfill) 
when A is smaller and/or when 7 is closer to t. 

Next, we proceed to solve the full lattice model numer- 
ically. On one hand, Fig. [1] displays the dressed quasi- 
particle spectrum, illustrating a band inversion at T > 0. 
On the other hand, Fig. [2] shows a topological crossover 
as a function of temperature. The topological "insula- 
tor" at T > is defined as a phase where m* < and 
\E* 



1 > 8T, the latter attribute being rather arbitrary. 
For many semiconductors, dE*/dT ~ 1 — 10 at high 
T [n| 13 24 1 ; hence it is conceivable that some real Dirac 
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FIG. 2: Phonon-induced topological quantum phase transi- 
tion (at T = 0) and topological crossover (at T > 0) for 
different values of the bare energy gap. The parameters 
(U, 7, A, t, u) c ) are the same as in Fig. [T] I: non-topological in- 
sulator (m* > 0); II: "crossover" regime (m* < and |m*| < 
4T); III: topological "insulator" (m* < and \m*\ > 4T). Be- 
cause ui c is fixed, the dashed lines in all the graphs converge 
to the same high-T asymptote (c.f. Eq. ©). 



insulators exist such that \E*\ > 8T deep in the inverted 
regime. In the "crossover" region, where surface states 
are emerging, m* < but \E*\ < 8T. In order to have 
a trivial insulator at T = turn into a topological "in- 
sulator" at high T, it helps if mju) c is larger because (i) 
at T = 0, U c increases for larger m (irrespective of ui c ) 
and (ii) at high T, \m*\ grows faster with T for smaller 
lo c (irrespective of m). 

Phonon-induced topological transitions in semimetals.- 
A simple lattice model for a 3D Dirac semimetal is 



h(k) = /i'°»(k) + A(7 2 l 2 



(8) 



where A > is an exchange field due to magnetic order 
that may originate e.g. from bulk dopants in a Dirac in- 
sulator. When m G (—A, A), the energy spectrum of h 
contains a pair of Weyl nodes of opposite chirality sep- 
arated in momentum space (along [001]): the system is 
a Weyl semimetal (WS) and has Fermi arcs [2j at the 
sample boundaries. For \m\ > A, the two Weyl nodes 
merge and the crystal turns into an insulator. Here we 
show that lattice vibrations, which preserve translational 
invariance only on average, can destroy or create Weyl 
nodes both at T = and T > 0. 

We evaluate E£ a by combining Eqs. (P, © and ©. 
The resulting phase diagram (Fig. [3]) evidences that, in 
presence of phonons, (i) A remains nearly unchanged, (ii) 
the WS phase is stable for m*(T) e (-A, A), and (iii) the 
Weyl nodes occur away from zero energy due to cf — > ^f- 
These generic observations are likely to hold for more 
elaborate models than Eq. ([8|). In addition, Fig. Oproves 
that a material which is insulating at T = may turn 
into WS at nonzero temperature. In order to observe this 
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FIG. 3: Phase diagram for Eq. ((8} in presence of electron- 
phonon interactions. The parameters (U, t, 7, A, uj c ) take the 
same values as in Fig. [T] with A = 20meV. For simplic- 
ity we have neglected any T-dependence of A by assuming 
that the magnetic order producing it is robust up to T ~ u c . 
When 771 > m>, increasing T may induce a WS phase; for 
m G (m<,m>), increasing T destroys the WS phase. When 
m < m<, the system is insulating for all T. The thick hori- 
zontal line segment is the window for the WS phase in absence 
of electron-phonon interactions and at T = 0. Insets: repre- 
sentative quasiparticle energy spectra for the insulating and 
WS phases, along [001]. Arrows indicate nodes with nonzero 
topological invariant. 



thermally-induced WS phase, the temperature at which 
the magnetic order giving rise to A dissolves must not be 
low compared to the Dcbye temperature. 

Conclusions- We have given a proof-of-principle for 
the impact of electron-phonon interactions on the band 
topology of insulators and semimetals, with an empha- 
sis on thermally driven topological crossovers. Although 
we have focused on phonons, other bosonic baths may 
likewise play a role in inducing or destroying topological 
phases in thermodynamic equilibrium. 

A necessary future task will consist of searching for 
real Dirac materials with appropriate band parameters 
to realize the qualitative effects discussed herein. Gener- 
ally, the more promising candidates will be soft crystals 
with strong electron-phonon coupling and a small zero- 
temperature bandgap that is nonetheless larger than the 
Debye temperature. Specifically, BiTl(Si_5 Se^ seems 
interesting because both the magnitude and sign of its 
bandgap are tunable However, it is not yet known 
how strong the phonon-induced bandgap rcnormalization 
is in this material. Should BiTl(Si_a Se^)2 become fer- 
romagnetic upon magnetic doping, there would also be 
an opportunity to alternate between insulating and Weyl 
semimetallic phases by tuning <5 and the temperature. 
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